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Abstract 

Dynamical systems driven by Gaussian noises have been considered 
extensively in modeling, simulation and theory. However, complex 
systems in engineering and science are often subject to non-Gaussian 
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fluctuations or uncertainties. A coupled dynamical system under non- 
Gaussian Levy noises is considered. After discussing cocycle prop- 
erty, stationary orbits and random attractors, a synchronization phe- 
nomenon is shown to occur, when the drift terms of the coupled system 
satisfy certain dissipativity and integrability conditions. The synchro- 
nization result implies that coupled dynamical systems share a dy- 
namical feature in some asymptotic sense. 

Key Words: Synchronization; Levy noise, Skorohod metric, ran- 
dom attractor, cadlag random dynamical system, impact of noise. 

Mathematics Subject Classifications (2000): 60H15, 60F10, 
60G17 

1 Introduction 

Synchronization of coupled dynamical systems is an unbiquitous phenomenon 
that has been observed in biology, physics and other areas. It concerns cou- 
pled dynamical systems that share a dynamical feature in an asymptotic 
sense. A descriptive account of its diversity of occurrence can be found in 
the recent book [33]. Recently Caraballo and Kloeden [6j [7] have proved 
that synchronization in coupled deterministic dissipative dynamical systems 
persists in the presence of various Gaussian noises (in terms of Brownian mo- 
tion), provided that appropriate concepts of random attractors and stochastic 
stationary solutions are used instead of their deterministic counterparts. 

In this paper we investigate a synchronization phenomenon for coupled 
dynamical systems driven by non-Gaussian noises (in terms of Levy motion). 
We show that couple dissipative systems exhibits synchronization for a class 
of Levy motions. 

Gaussian processes like Brownian motion have been widely used to model 
fluctuations in engineering and science. The sample paths of a particle driven 
by Brownian motion are continuous in time almost surely (i.e., no jumps), the 
mean square displacement increases linearly in time (i.e., normal diffusion), 
and the probability density function decays exponentially in space (i.e., light 
tail or exponential relaxation) [20]. But some complex phenomena involve 
non-Gaussian fluctuations, with peculiar properties such as anomalous dif- 
fusion (mean square displacement is a nonlinear power law of time) [I] and 
heavy tail (non-exponential relaxation) [35]. For instance, it has been argued 
that diffusion in a case of geophysical turbulence |28j is anomalous. Loosely 
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speaking, the diffusion process consists of a series of "pauses" , when the par- 
ticle is trapped by a coherent structure, and "flights" or "jumps" or other 
extreme events, when the particle moves in a jet flow. Moreover, anoma- 
lous electrical transport properties have been observed in some amorphous 
materials such as insulators, semiconductors and polymers, where transient 
current is asymptotically a power law function of time [2S1 d2] • Finally, some 
paleoclimatic data [H] indicates heavy tail distributions and some DNA data 
[28] shows long range power law decay for spatial correlation. 

Levy motions are thought to be appropriate models for non-Gaussian 
processes with jumps [25]. Let us recall that a Levy motion L(t), or L t , 
is a non-Gaussian process with independent and stationary increments, i.e., 
increments AL(t, At) = L(t + At) — L(t) are stationary (therefore AL has 
no statistical dependence on t) and independent for any non overlapping 
time lags At. Moreover, its sample paths are only continuous in probability, 
namely, P(|L(t) — L(t )\ > 5) — > as t — > to for any positive 5. With a 
suitable modification [1], these paths may be taken as cadlag, i.e., paths 
are continuous on the right and have limits on the left. This continuity is 
weaker than the usual continuity in time. In fact, a cadlag function has 
finite or at most countable discontinuities on any time interval (see, e.g., 
p. 118, [1]). This generalizes the Brownian motion Bit) or B t , since Bit) 
satisfies all these three conditions, but additionally, (i) almost every sample 
path of the Brownian motion is continuous in time in the usual sense, and 
(ii) the increments of Brownian motion are Gaussian distributed. 

This paper is organized as follows. We first recall some basic facts about 
stochastic differential equations (SDEs) driven by Levy noise in section 2, 
including a fact that the solution mappings of such SDEs generate random 
dynamical systems (RDS). In section 3, we formulate the problem of synchro- 
nization of stochastic dynamical systems driven by Levy noises. The main 
result (Theorem [T]) and an example are presented in section 4. 

2 Dynamical systems driven by Levy noises 

Dynamical systems driven by non-Gaussian Levy motions have attracted 
much attention recently [H[27]. Under certain conditions, the SDEs driven 
by Levy motion generate stochastic flows [H [18] , and also generate random 
dynamical systems (or cocycles) in the sense of Arnold (2]. Recently, exit 
time estimates have been investigated by Imkeller & Pavlyukevich [T3"l [T4] , 
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and Yang & Duan [53] for SDEs driven by Levy motion. This shows some 
qualitatively different dynamical behaviors between SDEs driven by Gaussian 
and non-Gaussian noises. 



2.1 Levy processes 

A Levy process or motion on M. d is characterized by a drift parameter 7 e M. d , 
a covariance d x d matrix A and a non-negative Borel measure v, defined on 
njd )) and concentrated on M d \ {0}, which satisfies 



/ 



(y 2 A 1) v(dy) < 00, 

Jm. d \{0) 

or equivalently 

y 2 

vidy) < 00. 



R d \{o> 1 + y 2 

This measure v is the so called the Levy jump measure of the Levy process 
L(t). Moreover Levy process L t has the following Levy-Ito decomposition 

L t = jt + B t + / xN(t, dx) + / xN(t, dx) (1) 

J\x\<l J\x\>l 

where N(dt, dx) is Poisson random measure, 

N(dt, dx) = N(dt, dx) - u(dx)dt (2) 

is the compensated Poisson random measure of L f , and B t is an independent 
Brownian motion on M. d with covariance matrix A (see [TJ [25], [2U E3])- We 
call (A, v, 7) the generating triplet. 

The next useful lemma provides states some important pathwise proper- 
ties of L t with two-sided time t 6 M. Here | • | denotes the usual Euclidean 
norm in M. d . 

Lemma 1. (Pathwise boundedness and convergence) 

Let L t be a two-sided Levy motion on M. d for which E|Li| < 00 and ELi = 7. 
Then 

(i) lim t - ±00 -L t = 7- 

6 
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(ii) The integrals e~ x ^~ s ^ dL s (uj) are pathwisely uniformly bounded in 
A > on finite time intervals [Ti,T2] in K. 

(Hi) The integrals f T e~ A ( <-s ) dL s (u) — > as A — > oo, pathwise on finite time 
intervals [T\,T2] in M. 

Proof, (i) This convergence result comes from the law of large numbers, in 
|25j . Theorem 36.5. 

(ii) Due to the continuous of function h(t) = e~ xt , on integrating by parts 
we obtain 

/ e~ xit - s) dL s (u) = L t {u) - A f e- x{t - s) L s (uj) ds. 

J ~oo J— oo 

Then we use (i) to conclude (ii). 

(iii) Integrating again by parts, it follows that 

f e- x ^ s) dL s (u) = (L t - L Tl )e- x ^ + X f e- x ^ s \L t (uj) - L s (u)) ds, 

from which the result follows. □ 

Remark 1. The assumptions on L t in the above lemma are satisfied by a wide 
class of Levy processes, for instance, the a-stable symmetric Levy motion on 
R d with 1 < a < 2. Indeed, in this case, we have ^ x]j>1 \x\v(dx) < oo, and 
then E\L\\ < oo; see f£W Theorem 25.3. 



Let us introduce the canonical sample space for Levy processes, the space 
Q = D(R,M. d ) of cadlag functions, i.e., continuous on the right and have 
limits on the left, defined on M. and taking values in M. d . 

If we use the usual compact-open metric, D(K, M. d ) is not separable. How- 
ever, it is complete and separable when endowed with the Skorohod metric 
[31 [2H], in which case we call D(M, lR d ) a Skorohod space. The Skorohod 
metric on £)(R, ~R d ) is defined as 

oo 

d ^ y)--=Y,^ 1A d °m( xm > y m )) f° r al1 x >v eD 

m=l 

where x m (t) := g m (t)x(t), y m (t) := g m (t)y(t) with 

1, if \t\ < m- 1 



9m(t) ■= < 



m — t, if m — 1 < |f| < m, 
0, if \t\ > m 
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and 

d m( x ,y) ■= \ni \ SU P 

AeA {-m<s<t<m 

where A denotes the set of strictly increasing, continuous functions from R 
to itself. 

Similarly, we can define a Skorohod space on a bounded time interval 
D([Ti, T2], R ). Then, in particular, the metric d° is the Skorohod metric on 
£>([-l,l],R d ). 

We recall the following compactness result (see [3], Page 116) in -D([7i, T2], M. d ). 

Lemma 2. (Ascoli-Arzela theorem in D([T\, T2], R rf ),) 
For 5 C [Ti,T 3 ], Zei 

w x (S) = sup{ \x(s) - x(t)\ : s,t e S} 
and for < 5 < 1 define 

w' x (5) = inf max w x ([ti-i,ti)) , 

ti 0<i<r 

where the infimum is taken over all the finite sets {t^} of points satisfying Ti 
— t < t\ < . . . < t r = T 2 with ti — ti-i < 5 for i — 1, 2 . . . , r. 

Then, a set B has compact closure in the Skorohod space D([Ti, T2], R d ) 
if and only if sup xeB sup t \x(t)\ < 00 and lim^o sn V x &A w 'x(^) = 0- 

2.2 SDE driven by Levy processes 

We consider the following stochastic differential equation (SDE) driven by 
Levy motion, which has continuous drift and Brownian motion components, 
namely 

dY{t) = b(Y(t-))dt + a(Y(t-))dBt + I F(Y(t—), x)N(dt, dx) 

J \x\<c 

+ f G(Y{t-),x)N{dt,dx) (3) 

J\x\>c 



log 



X(t) - X(s) 



V sup \x(t) - y(X(t))\ 



-m<t<m 
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where N(dt, dx) and N(dt, dx) are denned above, and the coefficients b, a, F, G 
are all assumed to be measurable. Here F and G may be different, while the 
positive parameter c may be different from 1, which allows greater generality. 
We introduce the d x d matrix 

a(x, y) = a(x)a(y f, x,y E R d , (4) 

and define 

d 

\\ a ( x ,y)\\ = X] \ a iA x ^y)\- 

i=l 

We make the following general assumptions for the SDE Q: 

A.l There exits K\ > such that, for all 2/1,1/2 G M. d , 

\KVi) ~ b (y2)\ 2 + \\a(yi,yi) - 2a(y u y 2 ) + a(y2,y 2 )\\ 
+ / \F{ yi ,x)-F{y 2 ,x)\ 2 u{dx) < K 1 \y 1 -y 2 \ 2 . 

J\x\<l 

A. 2 There exits K 2 > such that, for all y G M. d 

\b(y)\ 2 + \\a(y,y)\\+ [ \F(y, x)\\{dx) < K 2 \l + y\ 2 . 

J\x\<l 

A. 3 There exits 5 > 2 and K 3 > such that, for all yx, y 2 G M. d , 
[ \F( yi ,x) - F(y 2 ,x)\ p u(dx) < K 3 \ Vl - y 2 \ p , 

J\x\<l 

for all 2 < p < 5. 

From [Tj, Theorem 6.23, page 304, we have the following existence and 
uniqueness result for solutions of such SDE driven by Levy motion 

Lemma 3. Suppose that conditions (Al) and (A2) are satisfied and that 
the mapping y — > G(y,x) be continuous for all \x\ > c. Then there exists a 
unique global cadlag adapted solution of the SDE Q. 
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Note that a cadlag solution process has finite or at most countable dis- 
continuities on any time interval (see, e.g., p. 118, PQ). For more details 
about SDEs driven by Levy motions, see [HI [191 EB] • Due to the Levy-Ito 
decomposition ([I)), the following SDE, which we consider in the sequel, 

dY(t) = f(Y(t-))dt + g(Y(t-))dL t 

is a special case of ([3]). 

Remark 2. The reason to take the left limit in Y(t—) in the equation 
is to make sure that the cadlag solution process Y is predictable and unique 
I21f . For typographical convenience, however, we will write Y(t) instead of 
Y(t—) for the rest of this paper. Moreover, in the case of additive noise, i.e., 
if the noise intensity g(-) does not depend on the state Y, the distinction for 
left limit or not is not necessary, when we consider the integral form of the 
equation (jSJ), as f t f(Y(t—))dt = j to f(Y(t))dt for continuous f. In this 
case f(Y(t—)) has only countable discontinuous points and is thus Riemann 
and Lebesgue integrable. 

Remark 3. The above global assumptions do not hold for SDE, which we con- 
sider in the sequel, with a nonlinear dissipative drift f term such as x T f(x) 
< K — l\x\ 2 for some constants K > and I > 0. However analogous global 
existence and uniqueness results also hold in this case since the dissipativity 
condition prevents explosions and hence ensures otherwise local existence is 
global. See I2$j for more details. 

2.3 Random dynamical systems 

Following Arnold [2], a random dynamical system (RDS) on a probability 
space (Jl, P) consists of two ingredients: A driving flow 9 t on the prob- 
ability space Q, i.e., 9 t is a deterministic dynamical system; and a cocycle 
mapping <p : K x Q x Wt d — > M. d , namely, <p satisfies the conditions: 

<p(0, oS) = id R d, (p(t + s,u) — (pit, 6 s u) o (pis, to), 

for all uj G Q and all s,t G R. This cocycle is required to be at least 
measurable from the a-field B(R) x T x Bi^L d ) to the a-field Bi^L d ). 

For random dynamical systems driven by Levy noise we take Q = D(R, M. d ) 
with the Skorohod metric as the canonical sample space and denote by T := 
B(D(M., IR d )) the associated Borel cr-field. Let fii be the (Levy) probability 
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measure on T which is given by the distribution of a two-sided Levy process 
with paths in D(R, R d ). 

The driving system 9 = (8 t , t G R) on Q is defined by the shift 

(6 t uj)(s) := uj(t + s) - uj(t). (5) 

The map (t,u) —> 9 t u is continuous, thus measurable ([2] page 545), and the 
(Levy) probability measure is ^-invariant, i.e. 

^\A)) = fM L (A) 

for all A G J 7 , see p], page 325. 

Lemma 4. (RDS generated by SDEs driven by Levy motion) 
Suppose in addition to the assumptions of Lemma that condition (A3) is 
satisfied. Then there exists a unique cadldg adapted solution to and the 
solution mapping defines a RDS, which is continuous in x but cadldg in time. 

Proof. Let $ Sj t satisfy ([3]) with initial condition $ SjS (|/) = y, i.e. 

d$ Sit (y) = b($ 8 , t -(y))dt + a($ 8 , t -(y)))dB t + [ F($., t -(y)),a;)N(dt,dx) 

J\x\<c 

+ f G(Q 8 , t -(y)),x)N(dt t dx). (6) 

J\x\>c 

By Theorem 6.4.2 on page 322 and Corollary 6.4.11 on page 327 of pQ, $ is 
a Levy flow, and satisfies 

<5> 0tS+t {y,uj) = $o,t(&oM> d '(v))- 
We define tp : R x R d x Q -> R d by 

<p(t,y,w) = Q 0t t(y,u). (7) 

It follows that 

ip(t + s, y, w) = <p(t, ip(s, y, u), 8 (u)). 

Moreover, we note that <p(t,y,u) is continuous in y, measurable in u and 
cadlag in t. (cf. [I], page 336). It follows that the mapping (p is measurable 
from RxtfxOto R d . 

With this we only need Theorem 1.3.2 and Remark 1.3.3 in [2] (pages 
17-20) to complete the proof. □ 
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Remark 4. In view of Remark^ and analogous result holds for SDE with a 
nonlinear dissipative drift term fl^. Note that the perfection of crude discon- 
tinuous cocycles is considered in 13^ . 

We say that a family A = {A{uj)^uo G f2} of non-empty measurable 
compact subsets A(u) of M. d is invariant for a RDS (9,<p), if cp(t,cu, A(u>)) 
= A{9 t uo) for all t > and that it is a random attractor if in addition it is 
pathwise pullback attracting in the sense that 

H* d (ip(t, 9„ t uj, D(9- t uj)), A(uj)) -> as t -»• oo 

for all suitable families (called the attracting universe) of D = {D(u), uj G Q} 
of non-empty measurable bounded subsets D(oj) of M d , where is the 
Hausdorff semi-distance on M d . 

The following result about the existence of a random attractor may be 
proved similarly as in [301 El El ED EZ] • 

Lemma 5. (Random attractor for cadlag RDS) 

Let (0,(p) be an RDS on Q xl d and let (f be continuous in space, but cadlag 
in time. If there exits a family B = {B(u),uj G f2} of non-empty measurable 
compact subsets B{uj) ofW d and a > such that 

<p(t, 6- t u, D{B. t u)) c B(u), Vi > T b w , 

for all families D = {D(uj),u G Q} in a given attracting universe, then the 
RDS (6, ip) has a random attractor A = {A(u),u G Q} with the component 
subsets defined for each uj G Q by 



= f)\J<p(t,9- t u,B(9- t u)) 

s>0t>s 

Forevermore if the random attractor consists of singleton sets, i.e A(oS) = 
{X*(uu)} for some random variable X* , then X*(u) = X*{9 t uj) is a stationary 
stochastic process. 

We also need the following Gronwall's lemma from |22j . 

Lemma 6. Let x(t) satisfy the differential inequality 

^ x<g{t)x + h(t) 
at + 

where jfc.x := lj%|o+^f~^ right-hand derivative of x. Then 

x(t) < x(0)exp[ / g(r)dr] + / exp[ / g{r)dr]h{s)ds. 
Jo Jo J s 
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3 Dissipative synchronization 



Suppose we have two autonomous ordinary differential equations in M. d , 

where the vector fields / and g are sufficiently regular (e.g., different iable) 
to ensure the existence and uniqueness of local solutions, and additionally 
satisfy one-side dissipative Lipschitz conditions 

max {(xi - x 2 , /0i) - f{x 2 )), (xi - x 2 , g{x x ) - g{x 2 ))} < -l\x x - x 2 \ 2 

(9) 

on R for some I > 0. These dissipative Lipschitz conditions ensure existence 
and uniqueness of global solutions; see Remark [3] above. Each of the systems 
has a unique globally asymptotically stable equilibria, x and y, respectively 
|17j . Then, the coupled deterministic dynamical system in M. 2d 

^ = f(x)+\(y-x), ^L =g{x) + X (x-y) (10) 

with parameter A > also sastisfies a one-sided dissipative Lipschitz con- 
dition and, hence, also has a unique equilibrium (x x ,y x ), which is globally 
asymptotically stable [TTj . Moreover, (x x ,y x ) — > (~z,~z) as A — > oo, where 1 
is the unique globally asymptotically stable equilibrium of the "averaged" 
system in R d 

This phenomena is known as synchronization for the coupled deterministic 
system ( TTOi) . The parameter A often appears naturally in the context of the 
problem under consideraiton. For example in control theory it is a control 
parameter which can be chosen by the engineer, whereas in chemical reactions 
in thin layers separted by a membrane it is the reciprocal of the thickness of 
the layers, see [5]. 

Caraballo et al. (6j [7j showed that this synchronization phenomenon 
persists under Gaussian Brownian noise, provided that asymptotically stable 
stochastic stationary solutions are considered rather than asymptotically sta- 
ble steady state solutions. Recall that a stationary solution X* of a SDE sys- 
tem may be characterized as a stationary orbit of the corresponding random 
dynamical system (8, (/^(defined by the SDE system), namely, <p(t, uj, X*(cu)) 

= x*{e t u). 
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They considered a coupled system of stochastic differential equations 
(SDEs) in R 2d . 

dX t = [f(X t ) + X(Y t - X t )} dt + adBl 

(12) 

dY t = [g(Y t ) + \{X t - Y t )] dt + (3dB 2 . 

where a, (3 G R d are constant vectors with no components equal to zero, 
B], B 2 are independent two-sided scalar Brownian motions, and f,g satisfy 
the one-side dissipative Lipschitz conditions Q). This coupled system has 
a unique stationary solution (X t ,Y t ), which is pathwise globally asymp- 
totically stable. Moreover, the coupled system (Ti~2l is synchronized to the 
"averaged" SDE in R d 

dZ t = l - lf(Zt) + g(Z t )} dt + l -adB\ + ^(3dB 2 (13) 
in the sense that (X^, Y*) — > (Z™, Z°^) as A — > oo, where is the unique 



pathwise globally asymptotically stable stationary solution of ([13]) . 

The aim of this paper is to investigate synchronization under non-Gaussian 
Levy noise. In particular, we consider a coupled SDE system in W 1 , driven 
by non-Gaussian Levy noise in M. 2d 

U\X t ) + \{Y t -X t ))dt + adL\, 

= (g(Y t ) + X(X t -Y t ))dt + (3dL 2 t , 




(14) 



where a,/3,f,g are as above, and Lj,L 2 are independent two-sided scalar 
Levy processes satisfying conditions in Lemma [TJ We assume that this cou- 
pled system defines a random dynamical system ip (i.e.., it satisfies the as- 
sumptions in Lemma H] or some generalization of it). 

In addition to the one-side Lipschitz dissipative condition Q on the func- 
tions / and g, as in [B] we further assume the following integrability condition: 
There exists mo > such that for any m G (0,mo], and any cadlag function 
u : K — > M. d with sub-exponential growth it follows 

e ms \f(u{s))\ 2 ds < oo, / e ms \g(u{s))\ 2 ds < oo. (15) 

) J — oo 

Without loss of generality, we assume that the one-sided dissipative Lipschitz 
constant / < m . 
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In the next section we will show that the coupled system (I14p has a unique 
stationary solution (X t , Y t ) which is pathwise globally asymptotically sta- 
ble with (X t ,Y t ) — > (Z™,Z™) in Skorohod metric as A — > oo, pathwise 
on finite time-intervals [Ti,T 2 ], where Z t is the unique pathwise globally 
asymptotically stable stationary solution of the "averaged" SDE in M. d 

dZ t = l - [f(Z t )+g(Z t )} dt+ l -adL\ + l -(3dL 2 t . (16) 



4 Systems driven by Levy noise 



For the coupled system (TBj) . we have the follow two lemmas about its sta- 
tionary solutions. 

Lemma 7. (Existence of stationary solutions) 

If the Assumption ffl5l) holds, and f and g satisfy the one-side Lipschitz dis- 
sipative conditions 0, then the coupled stochastic system (|14|) has a unique 
stationary solution. 

Proof. First, the stationary solutions of the Langevin equations 

dX t = -XX t dt + a dLl dY t = -XY t dt + (3 dh\ (17) 

are given by 

X 4 A = ae~ xt I e Xs dL], Y t x = pe~ xt f e Xs dL 2 t (18) 

i-oo J— oo 

The differences of the solutions of (fl4|) and these stationary solutions satisfy a 
system of random ordinary differential equations, with right-hand derivative 
in time: 



^ (X t - X x ) = f(X t ) + \(Y t - X t ) + XX x , 
at + 

(19) 



j t + (Yt - Y x ) = g(Y t ) + X(X t - Y t ) + XY t } 



The equations ([19]) are equivalent to 
j t U x = f(X t )+X(V x -U x )+XY x ,j t V x = g(Y t )+X(U x -V x )+XX x (20) 
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where U x = X t - X x and V x = Y t - Y t x . Thus, 

\j t + (\ U t\ 2 + W t X \ 2 ) = {U x J(U x + X x )-f(X x )) + (V x ig (V x + Y x )-9(Y x )) 
+ (E/* f(X x ) + XY X ) + (V x , g(Y x ) + \X X ) - \\U X - V x \ 2 

< ~ W\ 2 + l^| 2 ) + j \f(X x ) + XY X \ 2 + \ \g(Y x ) + XX X \ 2 
Hence ,by Lemma [61 
\U t x \ 2 + \V x \ 2 < (\U x f + \V x f)e-^) 

+ ^ f<& [\f(X x ) + AF^| 2 + \g(Y x ) + \X X \ 2 ] ds 

This means that \U X \ 2 + \V X \ 2 is pathwise absorbed by the family B\ d = 
{B2 d (u),Lj G f2} of closed balls in M. 2d centred on the origin and of radius 
R\(lj), where the |-Ra(^)| 2 is given by 

1+_L1 / e ¥ [\f(X x (6 s uj)) + AF, A (^)| 2 + \g(Y x (6 s uj)) + \X x (6 s uj)\ 2 ] ds 

' J — oo 

Hence, by Lemma the coupled system has a random attractor A A = 
{A x (uj),uj G fi} with A a (cj) C B x d (u). 

However, the difference (AX t , AY t ) = {X} — X 2 , Y^ — Y 2 ) of any pair of 
solutions satisfies the system of random ordinary differential equations 





j- AX t 
dt + 


= fix]) 


- f(x 2 ) 


+ \(AY t -AX t ), 




I AY, 
dt+ 


= 9&?) 


~ 9(Y t 2 ) 


-X(AY t -AX t ), 


so 












f- A7 t 2 ) 


= 2(AX t: 


fW) - 


f(X 2 )) + 2(AY u g(Y t 1 )-g(Y 2 )) 



—2X\AX t — AY t \ 2 
< -2l(\AX t \ 2 + \AY t \ 2 ) 

from which we obtain 

\AX t \ 2 + \AY t \ 2 < (|AX | 2 + |AF | 2 ) e" 
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which means all solutions converge pathwise to each other as t — > oo. Thus 
the random attractor consists of a singleton set formed by an ordered pair of 
stationary processes (X t (a;), Y t (a;)). □ 

Remark 5. Using LemmaUl it can be shown that the random compact ab- 
sorbing balls B^J^oS) are contained in the common compact ball for A > 1. 

Lemma 8. (A property of stationary solutions) 

The stationary solutions of the coupled stochastic system (fT4l) have the fol- 
lowing asymptotic behavior: 



X X t (uj) - Y X t (uj) ^ as A 
pathwise on any bounded time-interval [Tx,T2] ofl 
Proof. Since 



d{X X t - Y X t ) 



-2X(X X -Y X ) + f(X X )-g(Y X : 



dt + adL\ — (3dL 



we have 



d(D£e 



X„2Xt\ 



2At 



f(X X )-g(Y X ] 



+ ae 2Xt dLl-pe 2Xt dL 



t ■ 



for with D x = X t — Y t , so pathwise 



< e - 2A(i " Tl) |D^| 



+ a 



-2A(t-s) 



Ti 



\f(K 



\g(Y x )\)ds 



f'e-W-) dL] 













By Lemma [T] we see that the radius R\(9 t u) is pathwise uniformly bounded 
on each bounded time-interval [Ti, T 2 ], so we see that the right hand of above 
inequality converge to as A — » oo pathwise on the bounded time-interval 
[Ti,T 2 ]. " □ 



We now present the main result of this paper. 
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Theorem 1. (Synchronization under non- Gaussian Levy noise) 
Suppose that the coupled stochastic system in ~R 2d 

dX t = (f(X t ) + X(Y t - X t )) dt + a dL\, 
dY t = (g(Y t ) + X(X t - Y t )) dt + (3dL 2 t , 



(21) 



defines a random dynamical system (9,ip). In addition, assume that f and 
g satisfy the integrability condition ffT5l) as well as the one-side Lipschitz 
dissipative condition ([9]). 

Then the coupled stochastic system (1211) is synchronized to a single averaged 
SDE in R d 



dZ t = \ [f(Z t ) + g(Z t )} dt + l -a dL] + ^ dL 2 



t ■ 



(22) 



in the sense that the stationary solutions of (1211) pathwise converge to that 
of ((22D, i.e. (X t ,Y t ) -> (Z?,Z?) in Skorohod metric on any bounded 
time-interval [Ti,T2] as parameter X — > oo. 

Proof. It is enough to demonstrate the result for any sequence A n — >• oo. 
Define 



t e 



Note that Z t A satisfies the equation 



dzl 



l - [f{X X t ) + g(Y")] dt + l -a dL] + ±(3 dL 



(23) 



(24) 



Also we define 



Z t := X t + Y tl te R, 
where X t and Y t are the (stationary) solutions of the Langevin equations 

dX t = -X t dt + adLl, dY t = -Y t dt + (3dL 2 t , (25) 

i.e. 



X t = ae~ 



! e s dL], Y t = pe- 1 f e s dL 



The difference Z^ — Z t satisfies pathwaise a random ordinary differential 
equation 



dt + [Zt 



t) 2 



(f{X X t )+g{Y X t ) 
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By Lemma [TJ we obtain 

\j t+ (Z*(u) -Z t (u))\ < ~\f(X*(u;))+g(Y*(u;))\ + l\X t { u ) + Y t (u>)\ 

< Mt^t^oS) < oo 

by the cadlag property of the solutions and the fact that these solutions 
belong to the common compact ball. We can use Lemma [2] to conclude that 
for any sequence A n — > oo, there is a random subsequence A nj (u;) — > oo, such 

that Z t ™ J (u) — Z t (u) — > Z£°(uj) — Z t {uj) in Skorohod metric as j — > oo. Thus 

Z t " J {oj) — > Z^°(uj) in the Skorohod metric as j — > oo. Now, 



M - ^ (<*>) = = ► 



( W ) - Z> ( W ) = it'M-^'M 



as A n . — > 00, so 



X> (c) = 2Z?> (u) - Y^ ( W ) - (u,), 
F> ( W ) = 2Z t ^ ( W ) - X> (u,) - Z t °», 



as A nj — > 00. Moreover, 



Z*-Z t = Z* 1 -Z Tl + ~£f(X X s )ds + ~£g(Y X s )ds 



1 ,* j ,t 



X s ds + - Y s ds 



2 .'ti 2 j Tl 



which converges pathwise to 
^° = Z% + ~£f(X^)ds+~j\(YT)ds 



1 /"^ . 1 



+z f - z Tl + - / x s ds + - / r s 



2 ./Ti 2 j Tl 



^ + s f f(K) ds+l f g(¥?) ds+^ f dh\ + t f dLl 
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on the interval [7i,T 2 ]. Therefore, is a solution of the averaged SDE 
for all t 6 1. The drift of this SDE satisfies the dissipative one-side 
condition (Q. It has a random attractor consisting of a singleton set formed 
by a stationary orbit, which must be equal to Z^°. 

Finally, we note that all possible subsequences of Z^ n have the same 
pathwise limit. Thus the full sequence Z A ™ converges to as A n — > oo. 
This completes the proof. □ 



4.1 An example 

Consider two scalar SDEs: 

dX t = -{X t + l)dt + dL], dY t = -{Y t + 3)dt + 2dL 

which we rewrite as 

dX t = -X t dt + dL\, dY t = -Y t dt + 2dLf, 

where L\ = 1 + L\ and L\ = 3/2 + L\. 

The corresponding coupled system (j2*Ij) is 

dX t = -X t dt + X(Y t - X t ) dt + dL\, 
dY t = -Y t dt + X(X t - Y t ) dt + 2dh\ 

with the stationary solutions 

-t rt 



X t x = I e- (A+1)( ^ s) coshA(t-s)^ + 2 f e- (A+1)(i - s) sinhA(t-s)^, 

J —oo J —oo 

Yj = f e- (A+1)( *- s) sinhA(t-s)^ + 2 f e- ix+1)(t ~ s) cosh X{t-s)dLl 

J —oo J —oo 

Let A — > oo, then 

where given by 

ZP°= [ le~ {t ~ s) dLl + f e~ {t - s) dLt 



2 



oo 



is the stationary solution of the following averaged SDE 

dZ t = -Z t dt+^ dh\ + dL\, 

18 



which is equivalent to the following SDE, in terms of the original Levy mo- 
tions L 1 and L 2 : 

dZ t = -{Z t + 2)dt+- dL] + dL 2 t . 

Acknowledgements. We would like to thank Peter Imkeller and Bjorn 
Schmalfuss for helpful discussions and comments. 
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